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, Abstract. We consider the problem of estimating the ensemble average of 

■ an observable on an ensemble of equally prepared identical quantum systems. 

We show that, among all kinds of measurements performed jointly on the 

q i copies, the optimal unbiased estimation is achieved by the usual procedure 

Othat consists in performing independent measurements of the observable on 
each system and averaging the measurement outcomes. 
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1. Introduction 

1 The astonishing precision of measurements currently available in quantum op- 

tics along with the growing demand of quantum devices of the new information 

■ technology |2 E] have revived the interest in the theory of quantum measurements 
Qh! 0] . The outcome statistics of a quantum measurement for all possible input states 
+!> 1 is described by a positive operator valued measure (POVM). The general optimiza- 
tion approach of Quantum Estimation Theory is to maximize over all possible 
POVM's an appropriate cost function, which depends on the context and on the 
specific use of the measurement. The output statistics can then be improved by 

I* ■ using multiple copies of the same quantum system, all prepared in the same state, 

. ,— i | and performing a suitable ensemble measurement over the copies. 

■ The experimental complexity of ensemble measurements is roughly classified by 
dividing them into three main categories: a) independent, b) separable, and c) en- 
tangled measurements. Category a) is described by tensor products of independent 
POVM's; b) by POVM's with separable elements only; c) by POVM's where some 
elements are entangled. Notice that the separability of POVM's generally does not 
correspond to a physical separability of measuring apparatuses 1 , and this classifi- 
cation remains essentially mathematical in nature. However, at least one can say 
that category b) contains all adaptive measurements (in which the choice of the 
measuring apparatus on the nth copy depends on the outcomes of previous mea- 
surements), whereas category c) contains those measurements that need quantum 
interactions between copies, implying that all copies during the measuring time 
must be at the same physical location, or, otherwise, that a "quantum memory" is 
available. 
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There exist separable measurements that cannot be performed by separate measuring appa- 
ratuses, i. ,e. by local operations and classical communication (LOCC) 
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Among the three categories of ensemble measurements, the category c) of entan- 
gled POVM's discloses the full exponential growth of the Hilbert space dimension 
versus the number of copies N for a virtually unlimited optimization of the statis- 
tical efficiency of the measurement, with the possibility of largely surpassing the 
performance of categories a) and b) [H| [7| E] ■ Indeed, over the last few years, it 
has been recognized that entangled measurements are usually more efficient than 
independent measurements, and the optimal measurement scheme is almost always 
entangled |1(JI II II fTSl . However, in some situations it has been also shown that 
asymptotically for N — > oo an equivalently optimal estimation may be achieved 
using just independent measurements over the copies [TH 1151 [TBI 1171 II 8| . 

In the above scenario it is natural to ask if the canonical procedure of averaging 
the outcomes of repeated measurements of an observable A over equally prepared 
systems is the best way of estimating the ensemble average (A) of A, or, instead, if 
a joint entangled measurement over the copies can improve the estimation. As we 
will see it turns out that the canonical procedure is indeed optimal, however, the 
derivation of this result is non trivial, and opens a general warning against easy 
assumptions and generalizations when evaluating statistical efficiencies of ensemble 
measurements. 

Let's be more precise, and fix precisely the scenario of the quantum estimation. 
Suppose one has a finite number N of equally prepared distinguishable identical d- 
dimcnsional quantum systems, which are described by the state p® N , and one wants 
to estimate the ensemble average (A) p = Tr[pA] of the observable A. Suppose now 
that one has unlimited technology at disposal, including measuring apparatus that 
can achieve any desired entangled POVM on all N systems jointly. The question is: 
which is the best measuring apparatus to choose in order to estimate {A) p with the 
minimum statistical error? What we will prove in the present paper is that the best 
estimation strategy is just the canonical procedure, which consists in averaging the 
outcomes of repeated measurements of the observable A over the equally prepared 
quantum systems. 

2. PERMUTATIONALLY INVARIANT POLARIZATION IDENTITIES 

In the derivation of our main result the following lemma will play a crucial role. 

Lemma 1. Any permutationally invariant operator X on H® N is completely de- 
termined by all ensemble averages Tr[X p® N ] on identical equally prepared systems. 

Proof. The statement of the lemma is equivalent to the following logical implica- 
tion 

(1) X EV N (H), \/peS{H), Tr[Xp® N } = => X = 0, 

where S(H) denotes the set of states on H, and 'Pat(H) the algebra of permuta- 
tionally invariant operators on H®^. Indeed, statement JTJ is equivalent to the 
statement that if Tr[Xp® N ] = Tr[Yp® N ] for all states p, then X = Y. 
Consider the following special states of the form 

N 

(2) Px = J2 AjI^XVj-I, a,- > 0, Xj ? Xi, i ? j, 

3=1 

with {ipj}N any set of N unequal states (not necessarily orthogonal). The trace 
Tr[Xpf N ] is a polynomial in Ylj=i Aj 3 , with YljLi x j — N and Xj > integers. 
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Now for Ti[Xpf N ] = all coefficients of the polynomial must vanish. In particular, 
the coefficient of YljLi Aj is given by 

(3) . . . (M^xul |Vi> • • • \^n) = o, 

where IT^ are the permutations of the N systems. By hypothesis we have U a XUl = 
X , then the vanishing of Ti[Xpf N ] for all states p\ implies 

(4) (^ 1 |...(Vjv|X|Vi)...|Vjv>=0, 

for all sets {ip^N- If we take \ipk) = a\<p) + p\<p±), by arbitraryness of a and j3 we 
have 

. . . <0| . . . <Viv|^|^i) ... |0> - - - |Viv> = 
{Vi|-..<^|.-.<Viv Wi>... |0±>... I ViV> = 

( } (-^il - - - <0^| - - - <Vjv|^r|-0i> ... |0> .. . |Vjv) = 

(Vi|...(^|...(Viv|X|Vi}...|</»i.)...|Viv) = 0. 

By repeating the same argument for different values of k and choosing <f) and (j)± 
as all possible elements of an orthonormal basis {4>j} we get 

(6) (<f> jl \...(ci> jN \X\ ( f )kl )...\<p kN )=0, V{ji},{*i}- 

Since all the matrix elements on an orthonormal basis of X are null, one has that 
X = 0. U 

Notice that the proof of the previous lemma contains the following interesting 
corollary 

Corollary 1. For any permutation invariant operator X on H® N the diagonal 
elements on factorized states completely determine X 

This is a kind of factorized polarization identity for permutation invariant oper- 
ators. 

3. The main result 

Let's now come back to the original problem of determining the optimal mea- 
surement for estimating the ensemble average of an observable. Consider a generic 
joint POVM P(r) on H®^, with outcome r providing an estimate of the expectation 
{A) p of the observable A on N identical systems all in the same state p. Clearly, 
one has A m ^ r ^ Am, with \ m and Am minimum and maximum eigenvalues of 
A, respectively. The POVM P(r) provides an estimate of the expectation (A) p if 
the conditional probability p{r\p) of estimating expectation value r for actual value 
Tr[Ap] is expressed via the Born rule as follows 

(7) p(r\p)dr = Tr[P(r)p® JV ]dr. 

Since the state p® N is permutation invariant, we can consider permutation invariant 
POVM's. Indeed, using invariance of p® N under permutations, one has 



P(r\ P ) dr=^E Tr[n ff p^4 P(r)] d r 



(8) 



Tr 



dr = Tr[n'(r)p ? 
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where the POVM 

(9) F 'M E ^E n ^ 

a 

is permutation invariant by construction. This means that for any POVM there 
is a permutation invariant one giving the same probability distributions for all 
states p® N . Therefore, without loss of generality, in the following we can assume 
that P{r) is permutation invariant. We will consider now the case in which the 
POVM is unbiased, namely that the averaging over r coincides with the value to be 
estimated. Mathematically this means that for all states p the following identity 
holds 

"Am 



(10) 



drrp(r\p) = Tr[Ap] 



The statistical error in the estimate is given by the r.m.s. of the probability distri- 
bution 



(11) 



e N (A) 



Am 



dr(r- (A) p ) 2 p(r\p) 



which for unbiased estimation equals 

"Am 



(12) 



e N (A) 



dr [r 2 p{r\p)] - (A) 



Since the only part which depends on the POVM is the conditional probability 
p{r\p), the optimization of the error resorts to minimize the quantity 

(13) J drr 2 Tr[P(r)p® N ], 

with the constraints 
(14) 



drP(r) 



(15) 

Using the following identity 



drrTr[P(r)p® N ] = (A) 



(16) 



(A) P = Tr 



(AT-l) 



V! 



)i4 



N ^ 



k=l 



with A^ = 7® (fe 1) (g) A (8) /® Af - fe ) by virtue of Lemma[I] we can recast Eq. {T5J) as 
follows 



(17) 



/ drrP(r) = -^A( fe )=9 



The operator A ^ defined as 
(18) A = 



drr 2 P(r) 
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allows to re-express the statistical error as follows 

(19) e N (A) 2 = Tr[Ap® N ]-(A) 2 p . 

In the representation in which A is diagonal, the constraints Ijl4(l and l|15(l become 

(20) jf " &rP{r)i k =5i k 

/■Am 

(21) / drrP(r)i k = eik 



whereas the error Ill2[l becomes 



Mm 

(22) e N (Af = J2(p® N ) nn / drr 2 P(r) ml - (A) 2 p . 

From Eqs. (|20l) and 1211) it follows that the diagonal elements P(r) nn are proba- 
bility densities versus r over [A m , Xm] , with average 8„„ , and upon denoting their 
variance by er 2 , we can write 

(23) e N {Af = Y,(P® N )nn(rt + OL) - {A%. 

n 

Therefore, €n(A) 2 is minimized by taking cr 2 = 0, corresponding to P(r) nn = 
i(9„„ — r). This implies that the outcomes of the optimal POVM are actually 
discrete, corresponding to r n = Q nn . In this discrete version, the POVM has 
P{r n ) nm = S nm [which also implies that 8„ m = 5 nm Q nn via Eq. I|21|l]. namely 
P{r n ) is projection- valued on the nth eigenvector of A, [when it happens that 
0„„ = mm for some m ^ n, then the projector has rank equal to the number of 
equal diagonal elements] . We have finally 

(24) e N (Af =J2(p® N )nnQL- {A? p . 

Tl 

Moreover, we have 

(25) Z = ^P( r „) 

n 

(26) e = ^e nn p(r n ) 

n 

(27) A = J2 & nn P ( r n) ■ 

n 

Since optimization makes 6 and A jointly diagonal, one has "}Z n Q 2 ln {p® N )nn = 
Tr[Q 2 p® N ], and using Eqs. I|ltj|) and l|17|l we can write the following expression for 
the minimal error 

(28) e N (Af = 1^ Tr[A®AWp» N ] - (A) 2 . 

Notice that the sum in the first term contains N terms with i = j equal to TrL4 2 p] 
and and N(N — 1) with i ^ j equal to (A) 2 , resulting in 



(29) e N (A) = 



N 
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namely the optimal error equals the statistical error occurring when measuring A 
separately on all the identical quantum systems in the state p, and then averaging. 
Indeed, the optimal POVM coincides with the spectral resolution of 8 = 4° ^2 n A^ n > 
on II v . 

Acknowledgements 

Giacomo Mauro D'Ariano and PP acknowledge financial support by the EC 
under the program ATESIT (Contract No. IST-2000-2968I), and by MIUR un- 
der programs Cofinanziamento 2003 and FIRB 2001. VG acknowledges financial 
support by MIUR through Cofinanziamento 2003. 

References 

[1] Quantum Optics, Experimental Gravity, and Measurement Theory, ed. P. Meystre and M. 
O. Scully, (Plenum Press, New York and London 1983) 

[2] Introduction to Quantum Computation and Information, ed. by H.-K. Lo, S. Popescu, and 
T. Spiller (World Scientific, Singapore, 1998). 

[3] I. L. Chuang and M. A. Nielsen, Quantum Information and Quantum Computation (Cam- 
bridge University Press, Cambridge, 2000). 

[4] P. Busch, P. J. Lahti, and P. Mittelstaedt, The Quantum Theory of Measurement, Lecture 
Notes in Physics Vol. 2 

[5] C. W. Helstrom, Quantum detection and estimation theory (Academic Press, New York, 
1976). 

[6] R. Derka, V. Buzek, and A. K. Ekert, Phys. Rev. Lett. 80, 1571 (1998). 
[7] A. Chefles, Phys. Rev. A 64, 062305 (2001). 

[8] G. M. D'Aria no, and P. Lo Prest i, M. G. A. Paris, Phys. Rev. Lett. 87, 270404 (2001). 
[9] A. S. Holevo, |quant-ph/0307225| (2003). 
[10] G. Chiribella, G. M. D'Ariano, P. Perinotti, and M. F. Sacchi, Phys. Rev. Lett. 93 180503 
(2004) 

[11] G. M. D'Ariano, C. Macchiavello, P. Perinotti, Phys. Rev. A (in press) 

[12] G. Chiribella, G. M. D'Ariano, P. Perinotti, and M. F. Sacchi, Phys. Rev A 70 062105 (2004) 

[13] G. M. D'Ariano, P. Perinotti, and M. F. Sacchi, Phys. Rev. A 72 042108 (2005) 

[14] E. Bagan, M. A. Ballester, R. Munoz-Tapia, O. Romero-Isart, Phys. Rev. Lett. 95 110504 

(2005) 

[15] M. A. Ballester, quant-ph/0506197 (2005) 

[16] E. Bagan, M. A. Ballester, R. Munoz-Tapia, O. Romero-Isart quant-ph/0505083 

[17] M. A. Ballester, Phys. Rev. A 70, 032310 ( 2004) 

[18] V. Giovannetti, S. Lloyd, and L. Maccone, quant-ph/0509179 (2005) 

1 QUIT Group, Dipartimento di Fisica "A. Volta", via Bassi 6, 27100 Pavia, Italy 

2 SCUOLA NORMALE SUPERIORE, PlAZZA DEI CAVALIERI 7, 56126 PlSA, ITALY 

a dariano@unipv. it, 'v.giovannettiOsns . it, c perinotti@f isicavolta.unipv. it, 



^http: / /www.qubit . it 



